Controlling the path of discretized light in waveguide lattices 
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A general method for flexible control of the path of discretized light beams in homogeneous waveguide 
lattices, based on longitudinal modulation of the coupling constant, is theoretically proposed. As 
compared to beam steering and refraction achievable in graded-index waveguide arrays, the proposed 
approach enables to synthesize rather arbitrary target paths. 
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Light propagation in waveguide lattices has received 
a great and continuous interest over the past few years 
QJII] , with the observation of a host of new phenomena 
such as optical Bloch oscillations [T] , Zener tunneling [5] , 
diffraction management [B] , dynamic localization [H] , 
Rabi oscillations [9 , Talbot imaging |10j . and Anderson 
localization just to mention a few. Lattice engi- 
neering enables to mold in a rather flexible way the flow 
of discretized light, hence providing altogether new op- 
portunities for applications [H IT2lfT4] . In spite of the 
discretized behavior imposed by the lattice, light trans- 
port in homogeneous waveguide lattices shear some com- 
mon features with optical beam propagation in homoge- 
neous media (hereafter referred to as continuous beam 
propagation). For instance, the path followed by a dis- 
cretized optical beam in a homogeneous waveguide lat- 
tice is straight like in a homogeneous medium, and beam 
spreading (diffraction) for both discretized and paraxial 
continuous beams is governed by the same algebraic law 
[T5] . A common method to control the path of continu- 
ous beams is to break the translational invariance of the 
medium in the direction transverse to the wave propaga- 
tion direction, realizing a graded-index (inhomogeneous) 
medium. The beam path is then determined by the pro- 
file of the refractive index according to the ray (eikonal) 
equation of Hamiltonian optics (see, for instance, [16]). 
A similar method holds for discretized beams. Graded- 
index waveguide lattices are usually realized by the in- 
troduction of an inhomogeneous profile of the propaga- 
tion constants or of the coupling strengths for the various 
waveguides HH E] , and can find applications in opti- 
cal steering and focusing [12,[T3]. For example, plasmonic 
aperiodic waveguide arrays have been recently proposed 
to realize deep sub- wavelength focusing and steering |14j . 
However, the synthesis of a graded-index structure that 
yields a desired beam path is a nontrivial issue. For beam 
propagation in continuous media, the ray (eikonal) equa- 
tion of Hamiltonian optics shows that beam steering and 
refraction can be realized by varying the refractive index 
along the longitudinal propagation direction, rather than 
in the transverse direction. This suggests that beam path 
control for discretized light beams might be realized in 
waveguide arrays without breaking the periodicity in the 
transverse direction. For example, steering of discrete 
optical solitons in optical lattices that fade away expo- 
nentially along the propagation direction was predicted 



in Ref. 18 . In this Brief Report we propose a simple and 
rather flexible method to control the path of a discretized 
beam in a homogeneous waveguide lattice, which is based 
on longitudinal modulation of the coupling constant. As 
compared to beam steering and refraction in inhomoge- 
neous graded-index arrays |14j . this approach enables to 
synthesize a rather arbitrary target path. 
Before discussing the beam path control method for dis- 
cretized light, let us recall the refraction properties of 
paraxial beams in a continuous medium with a refractive 
index n which varies solely along the paraxial direction 
z. For a monochromatic beam at wavelength A (in vac- 
uum), in the scalar approximation the complex electric 
field amplitude E(x, y, z) satisfies the Helmholtz equation 
d 2 E + V 2 E + k 2 n 2 (z)E — 0, where Vf is the transverse 
Laplacian and k — 2tt/X the wave number in vacuum. 
Assuming that n(z) varies slowly over one wavelength, 
by letting E(x,y,z) = tp(x,y, z) exp [ik f* d£n(£)], in 
the paraxial approximation the envelope ip satisfies the 
paraxial wave equation 



with a z-dependence of the diffraction strength. The so- 
lution ip(x, y, z) to Eq.(l), for an assigned initial field dis- 
tribution ip(x, y, 0) at the z — plane, can be simply ob- 
tained from the corresponding solution of Eq.(l) in vac- 
uum, in for n = 1. In fact, let 4>{x, y, z) be the solution to 
Eq.(l) with n = 1 and with 4>(x,y,0) = 2p(x,y 7 0). Then 
tp(x,y,z) — 4>(x,y,f Q d£/n(£)). Hence, since the cen- 
ter of mass of the beam in vacuum propagates along a 
straight path defined by the equations x(z) — x(0) + 8 x z, 
y(z) — y(0) + y z, where 9 X and 9 y are the paraxial 
beam angles at the z — reference plane, the center of 
mass for the beam ip(x,y,z) propagates along a curved 
path defined by equations x(z) — x(Q) + 9 X f Q d£/n(£), 
y(z) = y(0) + 6 y ci£/n(£). Thus a z-dependence of the 
diffraction strength for continuous beams results into a 
non-straight beam trajectory. Such a result can be ex- 
tended to discretized light. As the coupling rate n be- 
tween adjacent waveguides in a homogeneous array plays 
a similar role as the 'diffraction strength' for continuous 
beams, a longitudinal change of k is expected to curve 
the trajectory of the discretized beam. In fact, let us 
consider a rather standard tight-binding model describ- 
ing light transport at wavelength A in a homogeneous 
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waveguide array with lattice period a [TJ |2] and with a 
modulated (z- varying) coupling constant n(z) = Kof(z), 

ia n = -K f(z)(a n+1 + ct n -i), (2) 

where a„ is the modal amplitude of the light wave 
trapped in the n-th waveguide of the array, kq is a refer- 
ence value of the coupling constant, f(z) is a modulation 
function, and the dot stands for the derivative with re- 
spect to the longitudinal propagation distance z. The 
solution a n (z) to Eq.(2), for an assigned initial field dis- 
tribution a„(0) at the z = input plane, can be simply 
obtained from the corresponding solution of Eq.(2) for 
the non-modulated lattice, i.e. for f(z) = 1. In fact, 
let (j> n (z) be the solution to Eq.(2) with / = 1 and with 
0n(O) = On(0). Then a n [z) = <t> n (Jo d£f(0). Hence, 
since in a homogeneous array the center of mass of the 
beam <p n propagates along a straight path [15j . then the 
center of mass for the discretized beam a n {z) in the mod- 
ulated lattice propagates along a curved path. Indeed, 
assuming the normalization condition J2 n \ a n\ 2 = 1, af- 
ter letting (n)(z) = J2 n n \ a n(z)\ 2 for the beam center of 
mass, from Eq.(2) it readily follows that 

(n)(*) = (n)(0)-2/sop [' d£f(Q, (3) 
Jo 

where we have set p = Im a„(0)*a„_i(0)). The pa- 
rameter p entering in Eq.(3) is basically related to the 
tilting angle of the input beam. In fact, for a broad 
input beam tilted at the angle 9, one has a n (0) = 
|a„(0)| exp(ian) and thus p ~ — shier, where a = 
it9/8b and 9b = A/(2a) is the Bragg angle. From 
Eq.(3) it follows that a rather arbitrary path (n)(z) 
for the discretized beam can be achieved by a suitable 
choice of the modulation function f(z), namely f(z) = 
— (d(n)/dz)/(2K p). It should be noted that nonlinear 
propagation of discretized solitons in waveguide arrays 
with a longitudinally-modulated coupling constant was 
previously considered in Ref. [19] and shown to induce 
oscillations and decay of discrete solitons, however the 
possibility to exploit the longitudinal modulation to con- 
trol the path of beams in the linear propagation regime 
was not considered in such a previous work. Modulation 
of the coupling constant can be effectively realized by 
either waveguide axis bending or by out-of-phase mod- 
ulation of the propagation constants of adjacent waveg- 
uides (see, for instance, [7j I20H22] ). In such cases, the 
coupled-mode equations for the modal amplitudes c n of 
light trapped in the various waveguides read }21j 

i6 n = -K Q (c n+ i + Cn-x) + q n (z)c n , (4) 

where kq is the coupling constant between adjacent 
waveguides, and q n (z) — nG(z) for homogeneous arrays 
with axis bending, or q n (z) = (— l)"G(z)/2 for straight 
arrays with alternating modulation of the coupling con- 
stants. In the former case, the modulation function G(z) 
is related to the axis bending profile xq(z) by [7] 

G{z) — 2irn s Xo(z)a/ A, (5) 



where n s is the substrate refractive index, whereas in 
the latter case G(z) defines the alternating propaga- 
tion constant mismatch between adjacent waveguides. 
To establish an equivalence between the lattice mod- 
els (2) and (4), let us assume a sinusoidal modulation 
function G(z) with spatial period A and slowly-varying 
amplitude A, i.e. G{z) = A(z) cos(27rz/A), where A 
varies slowly over one spatial period A. After setting 
c„(z) = a„(z)exp[— i J* q n (Q], Eq.(4) can be cast in the 
equivalent form 

ia n = ~K n (z)a n+ i - < l _ 1 (z)a n _i, (6) 

where we have set n n (z) = Koexp[— i J* d£G(£)] 
in case of waveguide axis bending, or k„(z) = 
Koexp [(— l) n i f d£G(£)] in case of modulation of the 
propagation constants. Assuming that the spatial mod- 
ulation frequency f2 = 27r/A is larger than the coupling 
constant kq, at leading order in a perturbative analysis 
of Eqs.(6) [23] the evolution equations for the amplitudes 
a n take the form of Eqs.(2) with a modulation function 
f(z) given by 

where Jo is the Bessel function of first kind of zero or- 
der. In particular, if the envelope A(z) is varied such 
that A/Q, remains close to 2.405 (the first zero of Bessel 
Jo function), i.e. near the condition for suppression of 
evanescent tunneling [20], one has [24] 

f(z) ~ -0.52 - 2.405^ . (8) 

This means that the modulation envelope A(z) of axis 
bending or propagation constant is just mapped into the 
modulation of the coupling f(z), which is in turn related 
to the beam trajectory via the simple relation (3). There- 
fore, the synthesis of a quite arbitrary beam path can be 
realized in a very simple way. For instance, to realize a 
parabolic (cubic) path for a discretized beam, according 
to Eqs.(3) and (8) a linear (parabolic) change of A(z) is 
required. This is shown, as an example, in Figs. 1(a) and 
1(b). In the figures, the evolution of the discretized inten- 
sity distribution |c„(z)| 2 along the array, as obtained by a 
numerical analysis of Eqs.(4) for a tilted Gaussian input 
beam distribution and assuming a periodically-bent axis 
q n (z) = nA(z) cos(fiz), is depicted for a linear [Fig. 1(a)] 
and quadratic [Fig. 1(b)] variation of the amplitude A(z), 
resulting in an effective parabolic and cubic path of the 
beam, respectively. A similar result is obtained by con- 
sidering periodic out-of-phase modulation of the prop- 
agation constants rather than periodic waveguide axis 
bending. It is interesting to notice that the proposed 
scheme of beam path control also results in an effective 
engineering of the discrete diffraction. For single waveg- 
uide excitation at the input plane, the discrete diffraction 
pattern in a homogeneous lattice with a constant cou- 
pling Kg evolves according to |a„(z)| 2 = J^(2kqz), and 
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FIG. 1. (color online) Light evolution (snapshot of \c„\ 2 ver- 
sus the normalized propagation distance kqz) in a homoge- 
neous tight-binding lattice with a sinusoidally-bent axis with 
a normalized spatial frequency modulation fl/no — 5 and 
with amplitude A(z) = f2(2.405 — z/L) [in (a) and (c)], and 
A(z) = n(2A05-z 2 /L 2 ) [in (b) and (d)], where L is the array 
length. In (a) and (b) the array is excited by a tilted Gaus- 
sian beam c n (0) = exp(— n 2 /9) exp(i7rn/2), whereas (c) and 
(d) correspond to single waveguide excitation c n (0) = 5 n ,o- In 
(a) and (b) the beam follows a parabolic and cubic path, re- 
spectively. In (c) and (d) light spreading is faster than linear, 
corresponding to regimes of superdiffration. 
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FIG. 2. (color online) (a) Parabolic path of a broad Gaus- 
sian beam in a sinusoidally-curved homogeneous waveguide 
array with linearly- varying modulation amplitude. The fig- 
ure shows the evolution of the light intensity \ip(x, z)\ 2 versus 
propagation distance z for Gaussian input beam excitation 
(tilting angle 9 = 6b /2). The inset shows the behavior of 
the optical potential V(x) = n s —n{x) of the one-dimensional 
lattice. The values of other parameters are given in the text, 
(b) Same as (a), but for single waveguide excitation of the ar- 
ray. Note that light spreads quadratically versus propagation 
distance (superdiffractive regime), rather than linearly as in 
a homogeneous array with constant coupling. 



thus light spreads linearly with propagation distance fol- 
lowing the same ballistic transport law of electrons in 
tight-binding ordered crystals [25] . In the waveguide ar- 
ray with the modulated coupling constant, the spreading 
law can be engineered quite arbitrarily. For example, in 
case of linear or parabolic change of the amplitude A(z), 
regimes of superdiffraction can be realized, as shown in 
Figs. 1(c) and (d). The equivalence between the lattice 
models (2) and (4) has been established for a spatial mod- 
ulation frequency Q larger than the coupling constant kq, 
however it should be mentioned that beam steering can 
be achieved even for slow modulation frequencies. For 
the case of axis bending modulation, the expression of 
the beam path (n)(z) can be calculated in a closed form 
and reads (n) (z) = (n) (0) — 2n p d£f(£) + s(z), where 
f(z) is given by Eq.(7) and s(z) is a quasi-periodic func- 
tion with period A and with s(0) = 0. As for a fast 
modulation frequency s(z) — > according to the pertur- 
bative analysis, s(z) is non-negligible when f2 is of the 
order or smaller than kq. However, since at the planes 
z = A,2A, 3A, ... s(z) vanishes, a coarse beam steering 
control, at such discretized planes, can be realized even 



for slow spatial modulation frequencies. A fine beam 
path control requires, on the other hand, a fast modula- 
tion frequency. 

To check the feasibility of the beam steering method, nu- 
merical simulations of the full wave equation were per- 
formed. For the sake of definiteness, the case of axis 
bending modulation was considered. In the waveguide 
reference frame, the electric field envelope tf>(x, z) evolves 
according to the Schrodinger-type wave equation [3J [7] 

A 2 

iXd z Tp = -- — dlip + V(x)i> + n s x (z)xip, (9) 
2n s 

where A = A/(2?r) is the reduced wavelength of injected 
light, V(x) ~ n s — n(x) is the periodic lattice potential, 
n(x) is the refractive index profile of the array, and Xq(z) 
is the profile of axis bending. Equation (9) was integrated 
by a standard pseudospectral split-step method for a typ- 
ical waveguide lattice manufactured by femtosecond laser 
writing and excited at A = 633 nm [51 [52] . The refractive 
index profile of the array used in numerical simulations is 
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shown in the inset of Fig. 2(a). The corresponding band 
diagram, computed by a standard plane-wave expansion 
method, shows that the lowest band is well separated 
from the higher-order bands and its dispersion curve is 
well fitted by a sinusoidal curve. Hence, provided that 
the array is excited by an input beams tilted at an angle 
smaller than the Bragg angle and for a spatial modula- 
tion frequency small enough to avoid coupling to higher- 
order bands (radiation losses), the beam evolution in the 
lattice turns out be well described by the tight-binding 
lattice model (4). In the example shown in Fig. 2, the 
axis bending profile xq(z) has been chosen to realize a 
parabolic beam path as in Fig. 1(a), namely we assumed 
xq(z) = B(l — z/2L) cos(27rz/A), where L is the length of 
the waveguide array and B — 2.405 x AA/(47r 2 n s a) is the 
bending amplitude corresponding to tunneling inhibition 
(dynamic localization [7]). The spatial period A of the 
modulation used in the simulations is A = 3 mm, corre- 
sponding to B ~ 9 fiui for a bulk refractive index n s — 
1.42 and lattice period a = 9 /mi. Figure 2(a) shows the 
evolution of light intensity \ip{x, z)\ 2 along the array for a 
Gaussian-shaped input beam with spot size w = 21.6 /xm, 
tilted at half of the Bragg angle (0 = s /2 ~ 1°), i.e. 
ip(x,0) = exp(— x 2 /w 2 ) exp(iTTx/2a). According to the 
tight-binding model, the parabolic path followed by the 



beam is clearly visible, with negligible radiation losses in- 
duced by axis bending. Figure 2(b) shows the evolution 
of light intensity when a single waveguide of the array 
is excited in its fundamental mode at the input plane, 
leading to a superdiffraction regime for light spreading 
as in Fig. 1(c). 

To conclude, a flexible and simple method for the con- 
trol of the path of discretized light beams in homoge- 
neous waveguide arrays, based on longitudinal modu- 
lation of the coupling constant, has been theoretically 
proposed. As compared to beam steering and refraction 
control achievable in graded-index waveguide arrays, the 
proposed method enables to synthesize rather arbitrary 
target paths, and could be therefore of potential interest 
for beam steering applications in discrete photonics. Ow- 
ing to the quantum-optical analogy between light trans- 
port in waveguide arrays and coherent electronic or mat- 
ter wave transport in solid-state or matter wave systems 
[I] , the proposed method could be of interest beyond dis- 
crete optics. For example, it could be applied to control 
the path of coherent electronic wave packets as well as 
to realize superdiffusive coherent electronic transport in 
ac-driven quantum dot arrays |26j . 
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